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ABSTRACT 

As hundreds of gas giant planets have been discovered, we study how these planets form 
and evolve in different stellar environments, specifically in multiple stellar systems. In such 
systems, stellar companions may have a profound influence on gas giant planet formation 
and evolution via several dynamical effects such as truncation and perturbation. We select 
84 Kepler Objects of Interest (KOIs) with gas giant planet candidates. We obtain high- 
angular resolution images using telescopes with adaptive optics (AO) systems. Together 
with the AO data, we use archival radial velocity data and dynamical analysis to constrain 
the presence of stellar companions. We detect 59 stellar companions around 40 KOIs for 
which we develop methods of testing their physical association. These methods are based on 
color information and galactic stellar population statistics. We find evidence of suppressive 
planet formation within 20 AU by comparing stellar multiplicity. The stellar multiplicity 
rate for planet host stars is 0 (')q% within 20 AU. In comparison, the stellar multiplicity rate 
is 18%±2% for the control sample, i.e., field stars in the solar neighborhood. The stellar 
multiplicity rate for planet host stars is 34%±8% for separations between 20 and 200 AU, 
which is higher than the control sample at 12%±2%. Beyond 200 AU, stellar multiplicity 
rates are comparable between planet host stars and the control sample. We discuss the 
implications of the results to gas giant planet formation and evolution. 

Subject headings: 


1. INTRODUCTION 

Almost half of sun-like stars are members 
of bi nary or multip l e star systems (hereafter 
MSS , iRaghavan et al.l I20I0I; iDuauennov &: Mavod 
II99ID . During the formation of a star-planet 
system, the presence of a companion star is 
likely to perturb or truncate the protoplanetary 
disk (e.g., iHolman fc WiegertI 1 19991 I.Iang-Gondeill 
I2007D . These effects would affect the formation 
rate of planets around the stellar components of 
MSS. Indeed, observations of star forming regions 
show that protoplanetary disks around MSS have 
shorter lifetimes than disks around sing le star sys¬ 
tems (hereafter SSS. iKraus et al.l 1201211 . Eurther- 
more, after planet formation, dynamical interac¬ 
tions with stellar companions could affect the or¬ 
bital evolution of a planet and its survival rate 


via 1 

Derturbation (IWu & Murravl 120031 INaoz et al. 

2012 

), disk-driven migration (iLin & Paoaloizou 

1986 

: iTanaka et al.ll2002ll and nlanet-planet scatter- 

ing ( 

Rasio & Fordlll996l: IGhatteriee et al.ll2008l). 


Despite many theoretical and numerical studies of 
planets in MSS, the planet occurrence rate in MSS 
is still uncertain. There have been numerous es¬ 
timates of the planet occurrence rate for solar-type 
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stars based on the Kepler results (e.g.. lHoward et aH 
l20I2tlldessin et aDl20I3D : however, these studies do 
not distinguish between planet host stars in SSS 
or MSS. The lack of stellar multiplicity informa¬ 
tion prevents us from comparing the planet occur¬ 
rence rate between SSS and MSS. The comparison 
would provide insights into the influence of stel¬ 
lar companions on planet formation. The discovery 
of thous ands of planet candidates from the Kepler 
mission (iBorucki et al.l 120111 : iBatalha et al.l 1201,31 
iBurke et al.ll2014[l provides a unique opportuni ty to 
study planets in MSS (e.g.. lWang et al.ll2014aD 
There are two ways of studying planets in 
MSS. Eirst, one can select a sample of known 
MSS and then search for planets around each 
of the comp onent stars (the Planet-Quest ap¬ 
proach. e.g.. iKonackil l20^ lEggenberrer fc UdrvI 
l2007t iKonacki et al.ll2009HTovota et al.ll2009D . Um 
fortunately, for ground based surveys, the detec¬ 
tion efficiency of planets in MSS is affected by 
flux cont amination from the additional stellar com¬ 
ponents (jWright et al.l[20T^ . Flux contamination 
also reduces the s ignal to noise ratio (SNR) of a 
transit ing planet (iFressin et al.l 120131 iWang et al.l 
l2014al) . On the other hand, the Star-Quest ap¬ 
proach is an effective way to study planets in MSS. 
Using a sample of known stars with planets, it 
is possible to det ermine the fraction of MSS in 
that sample (e.g.. iLuhman fc ,Ia.yawardhaii^ 120021: 

Patience et al.l[^02t iWang et al.ll2014bl : iNgo et al.l 

20151) . From a technical standpoint, it is much eas¬ 
ier to detect stellar companions than planetary com¬ 
panions, making the Star-Quest approach much eas- 
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ier and more sensitive to lower mass planets. Pre¬ 
vious observational studies have suggested that the 
presence of a s tellar companion suppre sses planet 
formation le.g.. lEggenbererer et al.l[20TTI l. However, 
these studies have not fully considered biases in tar¬ 
get selection and planet detection or the observa¬ 
tional incompleteness in the se arch for stel l ar com- 
panions of planet host stars. iWang et al.l (l2014a[l 
discuss the above concerns and find that planet for¬ 
mation is suppressed by stellar companions with 
se parations up to 1500 AU. 

IWang et~al1 (j2014bl lah summarize the works fol¬ 
lowing the Star-Quest approach prior to 2014. Since 
then, more progress has been made. These new 
results suggest that visual stellar companions are 
not rare around pla net host stars. Using t he Lucky 
Imaging technique, iLillo-Box et al.l (l2014| l find that 
32.8% of Kepler planet host stars have at least 
one visual stellar companion within 6". Based on 
adaptive o ptics (AO) imag i ng of 87 Kepler planet 
host stars, [Dressing et al. (|2014ll find that 31.0% 
of planet host stars have stellar companions within 
4". iLaw et al.l (j2014ll observe 715 Kepler stars with 
planet candidates with the Robo-AO system. They 
find 7.4% of Kepler planet host stars have stel¬ 
lar companions within 2.5". However, the fraction 
of gravitationally bound companions is uncertain 
from t hese surveys. Using the speckle imaging tech¬ 
nique, iHorch et al.l (1201411 detect 49 stellar compan¬ 
ions within 1 arcsec around over 600 Kepler stars 
with planet candidates. The majority of the de¬ 
tected companions are likely to be gravitationally 
bound to the planet host stars based on a statisti¬ 
cal argument. Accounting for detection incomplete¬ 
ness, they conclude that the stellar multiplicity rate 
for planet host st ars is similar to star s in the so¬ 
lar neighborhood. iGilliland et al.l l)2014ll observe 23 
Kepler stars with small and cool planet candidates 
using the Hubble Space Telescope. They find evi¬ 
dence that physically-associated stellar companions 
are more common around planet host stars than 
around field stars in the solar neighborhood. The 
Kepler stars in the aforementioned studies have an 
average Kepler magnitude of 13.5. Assuming that 
they are solar-type stars, t he average distm rce is 
~600-700 pc. More recently. iNgo et all (j2015l l focus 
on a sample of host stars for hot Jupiters (HJs) that 
were discovered and confirmed by ground-based ob¬ 
servations. They use AO imaging and conduct 
multiple-epoch observations to confirm physical as¬ 
sociation by measuring common proper motions of 
stellar components. The stellar multiplicity rate for 
HJ host stars is almost twice as high as that for 
stars in the solar neighborhood for stellar separa¬ 
tions between 50 and 2000 AU. The overabundance 
of stellar companions for HJ host stars suggests the 
positive role of stellar companions in HJ formation 
and evolution. 

While rapid progress has been made since 2014, 
there are still a number of issues for the Star-Quest 
approach. First, there is a selection bias against 
MSS for work on planets detected by ground-based 
observations. MSS with small separations and con¬ 
siderable flux contamination are usually excluded in 


groun d-based surveys for planets (e.g.. lWright et al] 
l2012f) . The selection bias is difficult to quantify and 
correct, so converting the measured stellar multi- 
plici ty rate to planet o ccurrence rate is challeng¬ 
ing (|Wang et al.l 12014^ . Second, most studies fo¬ 
cus on a single detection technique for stellar com¬ 
panions. The high-angular resolution imaging tech¬ 
nique has been the dominant method. However, 
the physical association of stellar companions is dif¬ 
ficult to assess for Kepler stars. This becomes an 
issue when calculating the stellar multiplicity rate 
which concerns only gravitationally bound compan¬ 
ions. Third, the imaging techniques are not ef¬ 
fective in detecting stellar companions within or 
near the diffraction limits of telescopes, which cor¬ 
respond to 20-50 AU in physical separation. Stel¬ 
lar companions at smaller separations can be more 
effectively detected by measuring the radial veloc¬ 
ity (RV). So far, only a few studies following the 
Star-Quest approach use the RV technique in com¬ 
bination with the high-angular resolution imaging 
technique, which dramatically incre ases the search 
completeness for stellar companions (iKnutson et ahl 
IMa IWang et alJI2ni4al: iNgo et alJ[^15ll . 

Finally, the control sample to be compared is not 
a perfect sample. Field stars in the solar neigh¬ 
borhood usually serve as a control sample for stars 
without planets, but studies on planet occurrence 
rate suggest that the majority of sta r s host at least 


one 

planet fe.g., K 

avor et al.l 120111: IHoward et al. 

2012 

; iFressin et al.1 

20131: [Dressing & Charbonneau 

2013 

). In order to construct a more meaningful con- 


trol sample, we can find other stars that do not have 
planets down to a cert ain mass and up to a certai n 
orbital separation le.g.. lFggenberger fc Udrvll2007[l . 
or we can continue to use the field stars as a control 
sample but set a planet mass/radius and separation 
range to limit the level of contamination of planet 
host stars. For example, if we consider only gas gi¬ 
ant planets within ^1 AU, then the stars in the so¬ 
lar neighborhood can serve as a reasonable control 
sample because fewer than 10% of stars have gas 
giant planets with in 1 AU (iCumming et al.1 [20081 : 
iMavor et al.ll20li[l . 

To address the above issues, we have conducted 
a search for stellar companions for 84 Kepler stars 
with gas giant planets within ^1 AU. This sam¬ 
ple of stars is not biased against MSS because the 
Kepler mission doe s not apply a selec tion criterion 
that excludes MSS (|Brown et al.ll20lill . The typical 
full width at half maximum (FWHM) of images for 
Kepler target selection is 2.5" and thus these im¬ 
ages are ineffective at distinguishing binary stars. 
We use both RV data and high-angular resolution 
imaging data for these stars, so the search complete¬ 
ness is high compared to surveys employing only one 
technique. We assess the physical association of de¬ 
tected stellar companions using (1) their color in¬ 
formation and (2) a statistical argument based on a 
galactic stellar population model. We can compare 
the stellar multiplicity rate for this sample to that 
for field stars in the solar neighborhood without a 
significant contamination of planet host stars in the 
control sample. After these issues are resolved, we 
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can study the influence of stellar companions on gas 
giant planet formation with the sample 84 Kepler 
planet host stars. 

The paper is organized as follows. We describe 
the sample of Kepler stars with gas giant planets 
in ij2j AO observation and data reduction for these 
stars are presented in 321 We also discuss the physi¬ 
cal association of detected stellar companions in 321 
We synthesize the results of different techniques in 
the search for stellar companions in 321 These tech¬ 
niques include the RV and AO imaging techniques 
and the dynamical analysis. The stellar multiplicity 
rate of Kepler stars with gas giant planets is given 
in 33 Discussion and summary are given in 3Zl and 

33' 


2. SAMPLE DESCRIPTION 

From the NASA Exoplanet Archive^, we select 
Kepler Objects of Interest (KOIs) that satisfy the 
following criteria: (1), disposition of either Candi¬ 
date or Confirmed; (2), stellar effective temperature 
(Teff) lower than 6500 K; (3) stellar surface gravity 
(log^) higher than 4.0; (4), Kepler magnitude (Alp) 
brighter than 14th mag; (5), with at least one gas 
giant planet (3.8 < Rp < 22.0 i?©). In total, 

we select 84 KOIs with 97 gas giant planets. Stellar 
and orbital parameters for these KOIs are given in 
Table [2j The median distance of these KOIs is 580 
pc. There are 27 multi-planet systems among 84 
KOIs. 

There are 19 KOIs with RV observations. We ob¬ 
tain RV data from the Kepler Community Follow¬ 
up Observation Program® (CFOP). The major¬ 
ity of the RV data (14 out of 19 ) were taken 
with the HIRE S instrumen t (iVogt et al.l 119941) 
and reported in iMarcv et al.l (120 141) Exc eptions 
are KOI-1 (Tr ES^ lO’Don ovan et al.l 120061) . KOT 
3 (HA T-P-11 b.lBakos et al.ll2010l) . KOI-97 (Kepler- 
7 b. iLatham et all 120101) . KOI-128 (Kepler-15 
b. lEndl et al.l 2011 ). and KOI-135 (Kepler-43 b, 
iBonomo et al.l 2012). The Modified Julian Dates 
(MJDs) of the first and last RV data points and the 
number of RV data points for each KOI are given 
in Table [3 

Eor KOIs with high-angular resolution images 
from CFOP, we use the images to search for stellar 
companions. These images were taken at different 
telescopes including Keck, Palomar, MMT, Lick, 
and WIYN. For KOIs whose high-angular resolu¬ 
tion images are not available, we have taken AO im¬ 
ages using the PHARO (Paloma r High Angular Res¬ 
olution Observer) in strument (|Brandl et al.l Il997t 
iHavward et al.l[200111 at the Palo mar 200-inch tele¬ 
scope and the NIRC2 instrument (jWizinowich et aP 
[jniin i at the Keck H telescope. In total, we have 
taken AO images for 60 out of the 84 KOIs. Tele¬ 
scope and photometric band information is given in 
Table [3 The KOIs with AO images taken through 
this work are also indicated in Table [3 

3. AO OBSERVATION AND DATA REDUCTION 


^ http://exoplanetarchive.ipac.caltech.edu 
® https://cfop.ipac.caltech.edu 


3.1. AO Imaging with PHARO at Palomar 

We observed 40 KOIs i n the samp le with 
the PHARO inst rument (jBrandl et all Il997t 
IHavward et al.l I2001D at the Palomar 200-inch 
telescope. The observations were made between 
UT July 13rd and 17th in 2014 with seeing varying 
between 1.0" and 2.5". PHARO is behind the 
Palomar-3000 AO syste m, which provides a o n-sky 
Strehl of 86% in K band (jBurruss et al.l[20T^ . The 
pixel scale of PHARO is 25 mas pixel“^. With a 
mosaic IK xlK detector, the field of view (FOV) is 
25"x25". We normally obtained the first image in 
K band with a 5-point dither pattern, which had a 
throw of 2.5". The exposure time was set such that 
the peak flux of the KOI is at least 10,000 ADU for 
each frame, which is within the linear range of the 
detector. If a stellar companion was detected, we 
observed the KOI in J and H bands. 

3.2. AO Imaging with NIRC2 at Keck II 

We observed 27 KOIs in the sample w ith the 
NIRC2 instrument (jWizinowich et aH 120001 ) at the 
Keck H telescope. The observations were made on 
UT July 18th and August 18th in 2014 with excel¬ 
lent/good seeing between 0.3" to 0.8". NIRC2 is a 
near infrared imager designed for the Keck AO sys¬ 
tem. We selected the narrow camera mode, which 
has a pixel scale of 10 mas pixel“^. The FOV is 
thus 10"xl0" for a mosaic IK xlK detector. We 
started the observation in K band for each KOI. The 
exposure time setting is the same as the PHARO 
observation: we ensured that the peak flux is at 
least 10,000 ADU for each frame. We used a 3-point 
dither pattern with a throw of 2.5". We avoided the 
lower left quadrant in the dither pattern because it 
has a much higher instrumental noise than other 
3 quadrants on the detector. We continued obser¬ 
vations of a KOI in J and H bands if any stellar 
companions were found. 

3.3. Contrast Curve and Detections 

The raw data were processed using standard tech¬ 
niques to replace bad pixels, flat-field, subtract ther¬ 
mal background, align and co-add frames. We cal¬ 
culated the 5 -ct detection limit as follows. We de¬ 
fined a series of concentric annuli centering on the 
star. For the concentric annuli, we calculated the 
median and the standard deviation of flux for pix¬ 
els within these annuli. We used the value of five 
times the standard deviation above the median as 
the 5-cr detection limit. The median contrast curve 
and the 1-cr deviation of K band AO images we used 
in this paper are shown in Fig. [TJ Also plotted are 
detected stellar companions as indicated by aster¬ 
isks in Fig. [TJ These companions are brighter than 
the contrast curve, so the significance of detections 
is at least 5 a. In total, 59 stellar companions were 
detected around 40 KOIs. Their stellar and orbital 
properties are summarized in Tabled 

3.4. Comparison to Previous Work 

Among 59 stellar companions, 29 are newly de¬ 
tected around 22 KOIs in this study. Furthermore, 
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we add observations to 8 previously known stellar 
companions in additional color filters. The other 
stellar companions that were previously reported 
are noted with references in Table ID These obser¬ 
vation campaigns were carried out using a variety 
of instru ments at different telescopes, e.g.. AIRE S 
at MMT (|Adams et al.ll2012l: [Dressing et al.ll2014ll . 
PHARO at Palomar (I Adams et al.ll2012h . Robo-AO 
at Palor nar (iLaw et al.ll2014H . DSSI at WIYN and 
Gemi ni (iHorch et al.ll2014ll. and AstraLux at Calar 
Alto (jLillo-Box et alJl2014f l. When compared to 
previous work, we miss 11 stellar companions. They 
are marked with an asterisk in Table m All but 
one (KOI-372, AK = 4.0) stellar companions that 
we miss are very faint, with a differential magni¬ 
tude range between 7.2 and 8.2. Our pipeline does 
not identify these companions possibly because of 
different detection criteria. In one case (KOI-377, 
AJ = 6 . 8 ), the stellar companion is identified in K 
band, but not in J band. 


4. PHYSICAL ASSOCIATION 

For stellar companions detected by imaging tech¬ 
niques, we need to confirm that they are not op¬ 
tical doubles/multiples. Otherwise, the unassoci¬ 
ated stellar companions will systematically increase 
the stellar multiplicity rate and cause misinterpre¬ 
tations. To test physical association, the method 
of obtaining multiple-epoch images and measur¬ 
ing c ommon proper motion has been proven effec¬ 
tive (|Ngo et al.ll2015f) . In our case, Kepler stars are 
generally ^300-1000 pc away. While future obser¬ 
vations are scheduled, common proper motion mea¬ 
surements are relatively more difficult. Given only 
one epoch of observation, we can use color informa¬ 
tion of detected stellar companions and assess the 
proba bility of their physical association to pr imary 
stars (iLillo-Box et al.l[20l4lWang et al.ll2014all . De¬ 
tails of this approach are given in H4.ll For stel¬ 
lar companions with only single-band observations, 
color information is not available. We can assess the 
probability with a galactic stellar population simu¬ 
lation 1 H4.2D . 


4.1. Physical Association Based on Color 
Information 

We compare the distance of a KOI and its stellar 
companions. If their distances do not match within 
uncertainty, then they are likely to be optical dou¬ 
bles and the physical association is excluded. For 
the distan ce of a KOI , we fol low the method de¬ 
scribed in iWang et al.l (j2014a[ l. We calculate the 
distant modulus for each KOI. The V band appar¬ 
ent magnitude is obtained through the NASA Exo¬ 
planet Archive. The V band absolute magnitude is 
calculated u sing the Yale-Yonsei (Y 2) stellar evolu¬ 
tion model (jPemaraue et al.l 1200411 . The input pa¬ 
rameters for the Y2 model are Tgff, log 5 , age, and 
[Fe/H], which are also obtained through the NASA 
Exoplanet Archive. V band extinction (Ay) is ob¬ 
tained from the Mikulski Archive for Space Tele¬ 
scopes® (MAST). With the apparent and absolute 


http://archive.stsci.edu/ 


V band magnitudes and the extinction Ay, we can 
calculate the distance modulus of a KOI and thus 
its distance. For those KOIs whose extinctions are 
not available, we use AT band distance modulus as¬ 
suming zero extinction in K band. 

For stellar companions around a KOI, we use the 
color information, if available, to estimate their dis¬ 
tances. We have color information, i.e., multi-band 
detections, for 21 companions. We convert the dif¬ 
ferential magnitudes to the true color of the com¬ 
panion. Based on the color information, we estimate 
the effective temperature of a stellar comp anion us¬ 
ing Table 5 in iKraus &: Hillenbrandl (|2007li . For stel¬ 
lar companions detected in more than 2 bands, we 
use the mean effective temperatures weighted by un¬ 
certainties. Once the effective temperature is avail¬ 
able, we can find the corresponding K band absolute 
magnitude for a stellar companion. Its apparent AT 
band magnitude can be calculated from the differ¬ 
ential A' band magnitude and the apparent K band 
magnitude of the KOI. The AT band distance modu¬ 
lus can be calculated assuming zero extinction. The 
distance modulus can then be used to estimate the 
distance of a stellar companion. 

In the above calculation, extinctions in different 
bands need to be considered. Otherwise, a stellar 
companion would appear redder and closer. To ac¬ 
count for extinction in different band s, we use a lin¬ 
ear re lation between AxjAy and 1/A (jGordon et al.l 
I 2 OO 3 II . Since we are only interested in the wave¬ 
length region between 0.55 \im {V band) and 2.19 
/rm {K band), a linear relation is a reasonable ap¬ 
proximation. On one end, we assume K band ex¬ 
tinction to be zero. On the other end, we use the 
Ay from the MAST archive. Extinctions in r, i, z, 
J, and H bands are interpolated between Ay and 

Ak- 

The estimated distances of 21 companions are re¬ 
ported in Tabled For these companions with color 
information, 6 have estimated distances that are 2- 
cr inconsistent with the primary stars. Therefore, 
they are unlikely to be physically associated with 
the KOIs and thus are not considered in the follow¬ 
ing analyses. All 5 companions with color informa¬ 
tion and less than 1" angular separations have con¬ 
sistent distances with their KOIs. This is consistent 
with the finding that stellar companions with sub- 
arcsec se parations are most ly gravitationally bound 
to KOIs (|Horch et al.ll20l3l . For stellar companions 
with 1.0" to 3.0" separations, 2 out of 11 (18%) have 
inconsistent distances and are thus not physically 
associated with their KOIs. 

4.2. Physical Association Based on Galactic 
Stellar Population Model 

For the stellar companions without color infor¬ 
mation, we cannot adopt the method described in 
SH However, their physical association needs to 
be addressed because the frequency of optical dou¬ 
bles/multiples is not negligible: 6 out of 21 stellar 
companions with color information are not gravita¬ 
tionally bound. We therefore develop a statistical 
approach to assess the physical association of de¬ 
tected stellar companions. 
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Using the TRILEGAL galaxy 

model (|Girardi et al.l 1200511 . we run two sets 
of simulations. In the first set, we turn off binary 
parameters and calculate the fraction of optical 
doubles/multiples as a function of Ki, K 2 and A0, 
where Ki is the magnitude of primary star, K 2 is 
the magnitude of the brightest nearby star, A9 is 
the radius range in arcsec. In the second set of sim¬ 
ulation, we consider both optical doubles/multiples 
and gravitationally bound systems. From results 
of both sets of simulation, we can calculate the 
relative contribution of optical doubles/multiples 
and gravitationally bound stellar systems at a given 
combination of Ki, K 2 and A0, which allows us to 
calculate the probability of physical association in 
the absence of color information. 

In each simulation, ten fields with a FOV of 1 
square degree are simulated. These fields have dif¬ 
ferent galactic latitudes so the combination of the 
results from the fields gives a better statistical result 
of the entire Kepler FOV. We consider two different 
filters, J and K bands because all detections in sin¬ 
gle filter are in either J or K band. The majority 
(29 out of 38 ) are in K band. The physical associa¬ 
tion probabilities of detected stellar companions in 
single filter are given in Table U) We also provide a 
calculator for the probability of physical association 
as a function of Ki, K 2 and A9 in r, z, J, H and 
K filters^. 

4.3. Comparing Two Physical Association Methods 

We check the consistency of two methods for test¬ 
ing physical association. Since there are 21 stel¬ 
lar companions with color information, we can use 
this sample to perform the test: how physical as¬ 
sociation probabilities (in K band) correlate with 
acceptances/rejections based on color information. 
We divide the physical association probabilities into 
3 intervals, [0.00-0.33], [0.33-0.67], and [0.67,1.00]. 
For the lower probability interval [0.00-0.33], 2 out 
of 3 stellar companions (KOT98 and KOT377) are 
rejected based on color information at 2-a level. For 
the median probability interval [0.33-0.67], 2 out of 
3 stellar companions (KOT377 and KOT3444) are 
rejected based on color information. For the higher 
probability interval [0.67-1.00], 2 out of 15 stellar 
companions (KOT97 and KOT1812) are rejected 
based on color information. These results demon¬ 
strate consistency between these two methods. At 
a physical association probability smaller than 0.33, 
despite small number statistics, the majority (67%) 
of stellar companions are rejected based on color 
information. In contrast, at a physical association 
probability higher than 0.67, the majority of stel¬ 
lar companions (87%) show consistent colors to be 
physically associated. 

5. SYNTHESIZING AO OBSERVATIONS WITH OTHER 
TECHNIQUES 

While 59 stellar companions around 40 KOIs are 
detected via AO observations, the AO technique is 
not sensitive to stellar companions that are too close 

^ http://www.astro.yale.edu/jwang/CaLProb_PA.py 


to spatially resolve, nor is it sensitive to stellar com¬ 
panions that are too faint to detect with a sufficient 
SNR. By conducting simulations, we can calculate 
the search completeness of AO observations. 

We define a parameter space, a — i space, where 
a is the semi-major axis of a companion star, and 
i is the angle between the sky plane and the com¬ 
panion star orbital plane. We divide the parameter 
space into a grid (Aa = 0.5 AU, Ai = 10°). We 
simulate 1000 companion stars at each gridpoint in 
the a — i parameter space. The mass ratio distri¬ 
bution of simulat ed companions follows a Gau ssian 
distribution from iDuauennov &: Mavod (11991!) . i.e., 
q = m 2 /mi = 0.23, dq = 0.42. We use the median 
orbital eccentricity for binary stars (e = 0.4) and 
a random true anomaly distribution in simulations. 
If the contrast ratio (A Mag) between a simulated 
companion and the primary star is smaller than the 
value given by the 5 -ct AO contrast curve, then we 
record it as a detection. The median AO complete¬ 
ness contours are plotted in Fig. [2] 

For the parameter space on the a — i plane that 
AO is not sensitive to, we use other observations 
or techniques to constrain the presence of stellar 
companions. 

5.1. Radial Velocity Observation 

There are 19 KOIs in our sample with at least 3 
epochs of RV obser v ation. Following the descrip¬ 
tion of iWang et~aI1 (l2014a[) , we use the Keplerian 
Fittin g Made Easy (KFME) package ([Giguere et al.l 
l2012f) to analyze the RV data. For cases in which 
the number of RV data points are not adequate to 
constrain a Keplerian orbit, we use linear fitting 
to check if the RV data exhibit long-term trend. 
The RV data serve two purposes. First, they re¬ 
veal stellar companions via RV trends. Among 19 
KOIs with RV data, however, only KOT5 exhibits 
a RV trend. The stellar companion that can poten¬ 
tially i nduce the trend is c onstrained to be beyond 
7 AU (jWang et al.l 12014'^ . More recent RV data 
suggest that, in addition to two transiting planet 
candidates, two more distant components exist in 
KOT5 system (Howard Isaacson, private commu¬ 
nication). One is a sub-stellar companion with a 
period of ~2700 days and the other one is the AO- 
imaged stellar companion. Therefore, we consider 
the closest stellar companion to KOT5 to have a 
projected separation of 40.3 AU (Table Sj). 

The second purpose the RV data serve is to con¬ 
strain the presence of stellar companions in the non¬ 
detection cases. Given the RV data, we can study 
the completene ss of searching for s tellar companions 
by simulations (iWang et al.ll2ni4bO . Similar to the 
AO completeness study, we simulate 1000 compan¬ 
ion stars on each grid point and count the number 
of simulated companion stars that can be detected 
given the time baseline, observation epochs, and 
measurement uncertainties of the RV data. The me¬ 
dian RV completeness contours are plotted in Fig. 
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In addition to the RV and AO data, further 
constraints on potential stellar companions can 
be placed on multi-planet systems. There are 
27 (32% of the sample) multi-planet systems in 
our sample f or which we can apply a dynami¬ 
cal analysis (|Wang et al.l I2014bll . This dynam¬ 
ical analysis makes use of the co-planarity of 
multi -planet systems disco vered by the Kepler mis¬ 
sion (jLissauer et al.l 1201 If) . A stellar companion 
with high mutual inclination to the planetary or¬ 
bits would have perturbed the orbits and signifi¬ 
cantly reduced the co-planarity of planetary orbits, 
and hence the probability of multi-planet transits. 
Therefore, the fact that we have observed multiple 
transiting planet helps to exclude the possibility of 
a highly-inclined stellar companion. The dynamical 
analysis is complementary to the RV technique be¬ 
cause it is sensitive to stellar companions with large 
mutual inclinations to the planetary orbits. The 
parameter space to which the dynamical analysis is 
sensitive is shown in Fig. [2l 

5.3. Combining Results From Different Techniques 

For the RV and AO observations, detection com¬ 
pleteness contours are calculated based on simu¬ 
lations given the time baseline, cadence, measure¬ 
ment uncertainties, and the contrast curve. For 
the dynamical analysis, numerical integrations give 
the fraction of time when multiple planets can stay 
with small mutual inclinations (< 5°) s o that multi¬ 
ple tra nsiting planets can be observed (iWang et al.l 
l2014bli . We denote crVi cao and cda as the 
completenesses at a given point in the a — i pa¬ 
rameter space, overall completeness c is equal to 
1 — (1 — crv) X (1 — cao) X (1 — Cda)- 

The completeness is then integrated over the a — i 
parameter space. For the integration, distribution 
functions of a and i are necessary to account for 
contribution at different places in a — i parameter 
space. The result of the integration is sensitive to 
the adopted distribution function. Since the distri¬ 
bution function of a is uncertain for plant host stars 
and measuring the distribution is the main goal of 
this paper, we adopt an iterative approach to incor¬ 
porate a distribution of stellar companions. For the 
first it eration, we assume a log-no r mal distribution 
for a (|Duauennov fc MavodllQQit iRaghavan et al.l 
[MS). However, this distrib ution is not represe nta- 
tive for stars with planets (|Wang et al.ll2014all . so 
in the subsequent iterations we adopt the a distri¬ 
bution from )j6l The iteration stops when a distri¬ 
butions from two consecutive iterations differ less 
than 1% at any separations. 

We assume a random distribution of — cos i for 
systems with onl y one trans iting planet, and the i 
distribution from lHal^ (|1994f) for systems with mul¬ 
tiple transiting planets. The treatment for multiple 
transiti ng planet systems is detailed in iWang et al.l 
(j2014bl l. i.e., a coplanar distribution for stellar com¬ 
panions within 15 AU, a random — cosi distribution 
for stellar companions beyond 30 AU, and a mixture 
of the previous two i distributions for intermediate 
separations between 15 and 30 AU. 


5.4. Correcting For Detection Bias Against Planets 
in Multiple-Star Systems 

Planets in MSS are more difficult to find using 
the transit method because of flux contamination. 
The effect of this bias and a correction method have 
been discussed in IWang et~aI1 (l2014afl . We briefly 
introduce the method here. 

We conduct simulations to quantify the detec¬ 
tion bias against planets in MSS. For each KOI, we 
choose the one planet that gives the highest SNR. 
We add a companion star in the system and cal¬ 
culate the SNR in the presence of flux contamina¬ 
tion for two cases: planet transiting the primary 
star and planet transitin g the secondary star . If the 
SNR is higher than 7.1 ( Jenkins et al.l[2??T(il) . then 
the planet can still be detected, but with a lower 
significance. We randomly assign a stellar compan¬ 
ion (secondary star) to a KOI (primary star) and 
repeat this procedure 1000 times for both the pri¬ 
mary and the secondary star. We record the fraction 
of planet detections in 2000 simulations considering 
flux contamination. We designate the fraction to be 
a, which will be used in correcting for the bias of 
detecting planets in MSS. For example, a = 0.95 in¬ 
dicates that 95% of planets would still be detected 
in the presence of flux contamination. In order to 
account for the 5% missed planets, for every N MSS 
that host such a planet, we should use N /ato repre¬ 
sent the underlying MSS population that host such 
a planet. Since the transiting signal of gas giant 
planets is large, they are rarely missed in Kepler 
observations. Therefore, a is close to one in most 
cases. 

6. STELLAR MULTIPLICITY RATE FOR KEPLER 
STARS WITH GAS GIANT PLANETS 

The Kepler mission has provided us with a large 
sample of planet candidates. However, we do not 
know a priori whether a given planet host star is 
in SSS or MSS. Follow-up observations are critical 
in identifying additional stellar companions in plan¬ 
etary systems. Even in the case of non-detection 
with RV and AO, we can calculate the probability 
of a star being in a MSS based on the completeness 
study l T5.3p . For example, given the overall com¬ 
pleteness c and the stellar multiplicity rate (MR), 
the probability of the star having an undetected 
companion (or being in a MSS) within r (in AU) 
is: 

PM{a <r)= MR(a) / (1 — c{a,i))^{i) di da, 
J a<r J 

where uj(i) is a weighting function for i. For single 
planetary systems, ui{i)di = d{—cosi). For multi¬ 
ple planetary systems, uj(i)di is a piecewise function 
depending on stellar separation a ( 115.31) . The form 
of the weighting function for a, MR(a), was also 
discussed in 115.31 MR(a) is the a distribution of 
stellar companions for planet host stars. Here, we 
use MR(a) as a differential distribution, which is the 
derivative of a cumulative distribution MR(a < r), 
i.e., Equation|3l where both a and r are semi-major 
axis of an orbit. MR(a) and MR(a < r) are de- 
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rived in an iterative way. For each iteration, we use 
MR(a) from the previous iteration in Equation [T] 
to calculate pm , which is then fed into Equation [2] 
and [3] to calculate MR(a) in the new iteration. The 
iteration converges until MR(a) from the new iter¬ 
ation and MR(a) from the previous iteration agree 
within 1% at any separations. Following this proce¬ 
dure, we calculate the number of MSS, and the 
number of SSS, Ns- Since Nm and Ng are the sums 
of probabilities, they are not necessarily integers: 


N 

Nnia <r) = Y^{pM{a < r, k)/a{k)], Ns{a < r) = 


N 


Then a random number following the uniform dis¬ 
tribution between 0 and 1 is generated. If the ran¬ 
dom number is higher than the physical association 
probability, then the detection is excluded in the 
stellar multiplicity rate calculation. To account for 
the uncertainty in converting projected separation 
to semi-major axis, we assume randomly-orientated 
companion orbits. We use the median orbital eccen¬ 
tricity for companion stars (e = 0.4) and a random 
true anomaly distribution in simulations. The cal¬ 
culation for Nm, Ng, and the stellar multiplicity 
.rate is repeated for 1000 times for their values and 


( 2 ) 

where N is the total number of stars in the sam¬ 
ple, pM(k) is the probability of the fcth star being 
in a MSS, a{k) is the correction factor for the de¬ 
tection bias for planets in MSS (discussed in il5.4l) . 
Note that there is an implicit correction factor for 
SSS in Equation [21 but that this factor is 1. If a 
physically associated stellar companion is detected 
within a semi-major axis r to a KOI, then pM (a < r) 
is assigned to 1. We note that AO observation 
only measures projected separation. The conver¬ 
sion from projected separation to semi-major axis 
is addressed by a Monte-Carlo simulation assuming 
that stellar companions have randomly oriented or¬ 
bits (j ib.ID . We also assign a to 1 because no bias 
exists in this case since a planet has already been 
detected in a MSS. The cumulative stellar multi¬ 
plicity rate for planet host stars can be calculated: 


MR(a < r) 


Nuja < r) 

NM{a < r) -I- Ng{a < r) 


(3) 


6.1. Considering Physical Association Probability 
and Companion Orbital Orientation 

Not all AO detected stellar companions are phys¬ 
ically associated with the KOI. Therefore, we need 
to consider the probability of physical association 
when calculating Nm, which is later used for the cu¬ 
mulative stellar multiplicity rate calculation (Equa¬ 
tion [3]). Similarly, Nm may be different due to 
the orbit orientation of detected stellar compan¬ 
ions. AO observation only measures projected sepa¬ 
ration, but we need semi-major axis in Nm calcula¬ 
tion. Since orbital orientation, eccentricity and true 
anomaly are required to covert projected separation 
to semi-major axis, and these are not known for a 
single epoch AO observation, the conversion cannot 
be performed on an individual system. However, 
we can run a Monte-Carlo simulation to calculate 
Nm and its uncertainty due to physical association 
probability and companion orbital orientation. 

We developed two methods to calcu late the p rob- 
ability of physical association in 114.11 and 114.21 For 
detections in multiple filters, we estimate the dis¬ 
tance of a stellar companion based on its color in¬ 
formation. We exclude stellar companions whose 
distances are inconsistent with the KOI distance at 
more than 2-17 level. For detections in only one fil¬ 
ter, we estimate the probability of physical asso¬ 
ciation using a galactic stellar population model . 


6.2. Treatments For Different Stellar Separations 

For small separations, i.e., a < 10 AU, the RV 
data provide an effective constraint on stellar com¬ 
panions. As shown in Fig. [21 the completeness of 
the RV technique is higher than 50% for the ma¬ 
jority of parameter space within 10 AU. However, 
RV data are available for only 19 out of 84 KOIs. 
While considering all KOIs for stellar multiplicity 
rate within 10 AU seems to improve statistics, it 
in fact does not help because the majority of KOIs 
do not have data to constrain stellar companions 
within 10 AU. Instead, these KOIs without RV data 
outnumber KOIs with RV data and thus dominate 
the statistics. Since we use statistics of field stars in 
the solar neighborhood as an initial guess for stellar 
separation distribution for stellar companions (see 
ii5.3D . the lack of RV data for the majority of the 
KOI sample results in the lack of constraint for stel¬ 
lar companions within 10 AU. Therefore, the result¬ 
ing stellar separation distribution for these 84 KOIs 
would be similar to that of field stars in the solar 
neighborhood. The similarity of stellar separation 
distribution is not physical but rather a result of a 
lack of constraint from RV data for the majority of 
the sample. 

To avoid the above problem, we consider only 
KOIs with RV data when calculating stellar multi¬ 
plicity rate for small separations. To define small 
separations, we choose separations at which the 
completeness of AO data becomes higher than the 
completeness of RV data. Based on Fig. |21 the tran¬ 
sition separation is at 30-60 AU, so we adopt 50 AU 
as the transition separation. For stellar multiplicity 
rate within 50 AU, we consider 19 KOIs with both 
RV and AO data. For stellar multiplicity rate be¬ 
yond 50 AU, we consider all 84 KOIs for which we 
have AO data. 

One concern of using KOIs with RV data is the 
selection bias of RV observation and its potential 
influence on stellar multiplicity rate measurement. 
If RV observations are preferentially conducted for 
single stars or stars without significant flux contam¬ 
ination, then the stellar multiplicity rate for these 
KOIs would be lower because of selection bias. How- 
ev er, we have discuss ed this issue in Section 4.5 
of lWang et al.l (|2014a[) showing no evidence of such 
selection bias. For this work, we also checked the 
84 Kepler stars in our sample. For 19 stars that 
received RV follow-up observations, 5 have stellar 
companions within 2", one has severe flux contam- 
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ination, i.e., delta mag smaller than 2 mag. For 65 
stars without RV data, 11 have stellar companions 
within 2", 4 have severe flux contamination. The 
detection rates of stellar companions are compa¬ 
rable between stars receiving RV observations and 
stars without RV observations. Therefore, there is 
no evidence of selection bias of RV follow-up obser¬ 
vations, i.e., stars with RV data tend to have fewer 
stellar companions or fewer bright companions than 
stars without RV data. 

6.3. Stellar Multiplicity Rate vs. Stellar 
Companion Separation 

Fig. [3] shows the comparison between the 
cumulative stellar multipli city rate for field 
stars (blue hatched regio n, iDuauennov &: Mavod 
llQQHlRaghavan et al.ll2010ll and that for planet host 
stars (red hatched regions). The field stars serve 
as a control sample for comparison. Hatched re¬ 
gions represents l-cr uncertainties. For field stars 
in the so lar neighborhood, we adopt a 2 % un¬ 
certainty (|Raghavan et al.l [201011 . For planet host 
stars, we consider two sources of uncertainty. First, 
we conside r the uncertainty induced by physical as¬ 
sociation (ED- Second, we consider Poisson noise 
by propagating the uncertainty in Equation [H The 
two uncertainties are summed in quadrature for the 
final uncertainty. 

The stellar multiplicity rate for planet host stars 
is consistent with zero and stays flat at 01 q% for 
stellar separations smaller than 20 AU. For 20 AU 
< a < 200 AU, the stellar multiplicity rate is 
34%±8%. The increase of the cumulative stellar 
multiplicity rate from [0-20] AU to [20-200] AU is 
much faster than the control sample whose cumu¬ 
lative stellar multiplicity rate changes from 18% to 
30%. Beyond 200 AU, the stellar multiplicity rate 
for planet host stars is higher than that for the con¬ 
trol sample, but they are consistent within measure¬ 
ment uncertainty. The slopes for the cumulative 
stellar multiplicity rate change are similar between 
planet host stars and the control sample. 

7. DISCUSSION 

7.1. Interpretation of the Stellar Multiplicity of 
Stars With Gas Giant Planets 

According to Fig. [3l there may be three stellar 
separation ranges in which stellar companions affect 
gas giant giants formation and evolution differently. 
For stellar separations smaller than 20 AU, planet 
formation is suppressed. No stellar companion has 
been found within 20 AU for Kepler stars with gas 
giant planets. This leads to a zero stellar multiplic¬ 
ity ( 0 ^ 0 %) that is significantly lower than that for 
the control sample (18%±2%). 

For separations between 20 AU and 200 AU, we 
notice a drastic increase of the cumulative stellar 
multiplicity rate for planet host stars. In contrast 
to the stellar multiplicity rate of 12 %± 2 % in this 
separation range for the control sample, the stellar 
multiplicity rate is 34%±8% for planet host stars. 
The higher stellar multiplicity rate for planet host 
stars suggests that stellar companions in this sep¬ 


aration range play an important role in gas giant 
planet migration, e.g., via the Lidov-Kozai mecha¬ 
nism. Therefore, for 20 AU < a < 200 AU, the role 
of stellar companions in planet migration is more 
important than their role in suppressing planet for¬ 
mation. However, there may be a caveat assert¬ 
ing the role of stellar companions in planet migra¬ 
tion. The search completeness is ^^30-50% between 
20 and 200 AU (Fig. [2|), which is the lowest in the 
a — i parameter space. Therefore, the stellar mul¬ 
tiplicity rate in this separation range is the most 
uncertain. 

For separation beyond 200 AU, the stellar mul¬ 
tiplicity rate for planet host stars is comparable to 
that for the control sample, which indicates that gas 
giant planet formation and evolution is not signifi¬ 
cantly affected by a stellar companion. 

7.2. Comparison to Stars with Small Planets 

iWang et al.l (|2014all measured the cumulative 
stellar multiplicity rate for a sample of planet host 
stars. This sample is dominated by stars with 
smaller planets, 43 out 56 stars have only transit¬ 
ing planets that are smaller than 3.8 R m. We com¬ 
pare th e stellar multiplicity rates from IWang et al.l 
(|2014all and this work, which focuses on gas gi¬ 
ant planet host stars. While they are qualitatively 
identical, the characteristic separations are differ¬ 
ent. For example, the stellar multiplicity rate for 
small planet host stars intersects with control sam¬ 
ple at ^1500 AU whereas the intersection takes 
place at ^100 AU for gas giant planets. The ef¬ 
fective range of a stellar companion is an order of 
magnitude larger for small planets than for large 
planets: smaller planets are more prone to the in¬ 
fluence of a stellar companion. 

There may be several explanations. First, in 
a multi-planet system, the timescale for pericen- 
ter precession and nodal precession increase s with 
decreasing planet mass (jTakeda et al.l [200^ . For 
planet systems of the same orbital configurations, 
the Kozai timescale is more likely to be shorter than 
the timescale for pericenter precession and nodal 
precession for smaller planets, which makes these 
systems dominated by the Kozai effect due to the 
stellar companion. Therefore, smaller planets are 
more prone to the influence of a distant stellar com¬ 
panion because of a weaker planet-planet dynami¬ 
cal coupling. Second, for planet systems with both 
small and large planets, dynamical interaction be¬ 
tween these two types of planets tend to eject more 
small planets than large planets (Xie et al. 2015 in 
prep.). In the presence of stellar companions, the 
dynamical interaction becomes more frequent and 
thus leads to higher loss of small planets. 

7.3. Correlation Between Stellar Companion 

Properties and Planet Properties 

Planet formation is subject to the influence of stel¬ 
lar companions. Therefore, planets with different 
properties (e.g., orbital period and planet radius) 
may be a result of different properties of stellar com¬ 
panions. Here, we study whether the difference of 
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stellar companions results in the difference of planet 
properties. 

For properties of stellar companions, we focus on 
their differential magnitudes (A Mag). Since the 
majority (51 out of 59 ) of detected stellar compan¬ 
ions have differential magnitudes in K band, we use 
the differential magnitudes in K band. For those 
whose K band differential magnitudes are not avail¬ 
able, we do not use them in the following analysis. 

To find evidence that stellar companion properties 
affect planet properties, we adopt a K-S-test-based 
method that has been used to search for different 
popu l ations divided by a p arameter (|Quinn et al.l 
120141 iBuchhave et al.l 120141) . The method is de¬ 
scribed as follows. We choose a value x for a planet 
property parameter. The value divides the sample 
into two sub-samples. We then compare the stellar 
companion properties of two sub-samples with the 
K-S test. The p value of the K-S test is recorded. 
By varying x and repeating the K-S test, we ob¬ 
tain a function p{x). If there are certain x val¬ 
ues that result in significant difference between two 
sub-samples, these values may represent character¬ 
istic values that divide different planet populations. 
These populations are a result of different proper¬ 
ties of stellar companions. Because not all stellar 
companions are physically associated, we need to 
account for the effect of inclusion of optical dou¬ 
ble/multiples. The uncertainty due to this effect 
is addressed in a Monto Carlo simulation. In each 
trial, we draw a subset of detected stellar compan¬ 
ions based on the probability of their physical as¬ 
sociation. We apply the K-S-test-based method to 
the subset and record p{x). We repeat the process 
for 1000 times and find the median and confidence 
intervals at different levels. 

Fig. U shows the p value in K-S test as a function 
of planet orbital period. The K-S test compares 
the A Mag of two sub-samples divided by orbital 
period. The K-S test-based method suggests that 
there may be three different populations of gas gi¬ 
ant planets. We caution that this finding is not 
conclusive at this stage given small number statis¬ 
tics. If using p=0.05 as a threshold, the dividing 
orbital periods are ~10 and ~70 days. The p-value 
dip at P^IO day is broad extending from 7 to 22 
days. Hereafter we use P~10 days as a represen¬ 
tative value. The stellar properties of these three 
populations may be distinctively different. Fig. [5] 
shows the A Mag distribution. Stars with gas giant 
planets with P < 10 days tend to have more stellar 
companions with small differential magnitudes (A 
Mag < 2). In comparison, A Mag for stars with gas 
giant planets with P > 70 days tend to be fainter 
but peak at A Mag ^ 3. However, the peak may be 
due to a lack of sensitivity for fainter stellar compan¬ 
ions. The A Mag for Stars with gas giant planets 
with 10 < P < 70 days lies in between the previous 
two populations. 

We check whether the finding of three distinct 
populations is prone to inclusion of false positives. 
One major cause of false positive is und erestimation 
of rad ius due to flux contamination. iHorch et al.l 
(|2014f) estimate the effect as a function of A Mag 


in Kepler band in two scenarios: object transiting 
primary star and secondary star. After applying 
their calculation, there are five KOIs that could 
be false positives. All the potential false positives 
are for the scenarios in which an object is transit¬ 
ing secondary star. These KOIs are KOI-17, KOI- 
375, KOI-633, KOI-2672, and KOI-3678. Among 
them, KO I-17, also known as Kepler-6, has been 
confirmed (|Dunham et al.l[201Clf) . KOI-2672 is a sys¬ 
tem with multiple transiting planet candidates, so 
the f alse positive probability is low (jLissauer et al.l 
l2012f l. The stellar companions for KOI-375, KOI- 
633, and KOI-3678 are only detected in K band 
and have differential magnitudes of 3.3, 3.9 and 3.3 
mag. If physically associated, the K band differen¬ 
tial magnitudes would put the stellar companions 
in the range of early-type M dwarfs. Such low- 
mass stars are less likely to harbor gas giant planets, 
which makes the scenario in which an object tran¬ 
sits the secondary star le ss probable. Given t he false 
positive rate of ^^20% (iFressin et al.l 1201^ . there 
would be at most one KOI that is a false positive, 
which will not significantly reduce the peak seen in 
Fig. [5l However, we must caution again that the 
three-population hypothesis is preliminary and re¬ 
mains to stand the test of more observations and a 
larger sample. 

7.4. Comparison to Other Results 

The stellar multiplicity rate of planet host stars 
has been a research interest for a number of works 
(see 31] for references). All previous results concern¬ 
ing the Kepler sample were for small planet host 
stars. This is because Kepler discoveries are domi¬ 
nated by small planets and these planets are of great 
intere st for potential habitable worlds. iNgo et id] 
(|2015f) used a sample of stars with gas giant planets 
that were discovered by ground-based RV and tran¬ 
siting surveys. In their “Friends of Hot Jupiters” 
project, they found that the stellar multiplicity for 

50 stars with hot Jupiters (HJs) is 48 ± 9% for 
the stellar separation range of 50-2000 AU. We 
have 28 HJs in our sample, the stellar multiplic¬ 
ity rate we find for these separations is 25 ± 20%. 
However, if considering separations smaller than 50 
AU, the stellar multiplicity rate within 2000 AU is 

51 ± 13%. The significant difference is due to the 
drastic change of stellar multiplicity rate between 
20 and 200 AU (Fig. 0. 

In 47.31 we suggest that the stellar companions 
around HJ host stars are on average brighter than 
stellar companions around stars hosting longer- 
period gas giant planets (Fig. [5]). In our sam¬ 
ple, all stellar companions around HJ host stars 
have K band differential magnitudes smaller than 
2.5 with one exception of KOI- 17 (A Ma. g = 3.8 
mag). The trend is absent in iNgo et aP (1201511 . 
None of their detected stellar companions have dif¬ 
ferential magnitudes smaller than 2.5 mag. One 
explanation is that they focus on stars with HJs 
discovered by ground-based RV and transiting sur¬ 
veys. In these surveys, a strong bias against MSS 
exists in the target sele ction and f ollow- up obser¬ 
vations. The fact that INgo et ahl (|2015ll detected 
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many fainter stellar companions suggests that there 
may be a missing population of HJ systems in equal- 
mass binary systems. Recent discovery of HJs in 
WAS P-94 AB system is one exam ple of such sys- 
te ms (iNeveu-VanMalle et'mi2014ll . 

iLaw et al.l ()2014ll observed 715 KOls using the 
Robo-AO system, they found that stars hosting 
short-period (P < 15 days) giant planets are 2-3 
times more likely to have stellar companions than 
their longer-period counterparts. For our sample, 
we have 32 giant planets with period shorter than 
15 days. When comparing the stellar multiplicity 
rate for stars hosting short-period planets and the 
rest of our sample, we do not find a signihcant dif¬ 
ference. In fact, the stellar multiplicity rates for the 
stars with short-period planets, stars with planets 
with P > 15 days, and the entire sample are all 
^50% within 5000 AU. For stars with short-period 
planets, they tend to have brighter stellar compan¬ 
ions, which are easier for Robo-AO to detect. For 
stars with planets with longer periods, their stel¬ 
lar companions tend to be fainter with differential 
magnitudes peaking at A Mag = 3.0 in K band. 
These faint stellar companions are 1-2 mag fainter 
in the visible bands at which Robo-AO operates. 
Thus they are likely to be missed by Robo-AO given 
the detection limits are ^4-5 mag in r and i band 
for the median performance. The portion of fainter 
stellar companions missing in the Robo-AO survey 
may explain their finding that stars hosting short- 
period (P < 15 days) giant planets are 2-3 times 
more likely to have stellar c ompanions. Ther efore, 
we suspect that the finding in lLaw et al.l (|2014ll may 
be caused by a lack of sensitivity to faint stellar 
co mpanions. 

iJang-Condelll (|2015f) quantified planet formation 
efficiency in close binaries. The paper gave an em¬ 
pirical equation of gas giant planet formation ef¬ 
ficiency as a function of primary star mass, mass 
ratio, binary separation and binary orbital eccen¬ 
tricity. Our finding of a suppressive planet forma¬ 
tion within_20_AU is con sistent with t he conclusion 
in iJang-Condelll (|2015ll . However, iJang-Condelll 
(|2015ll predicted that gas giant planets in close bi¬ 
naries with separations smaller than 20 AU are not 
entirely impossible. In comparison, there is no such 
planetary system in our sample. The closest stel¬ 
lar companion we have detected is at 40.3 AU pro¬ 
jected separation. There are several other plane¬ 
tary systems in close binaries detected by the RV 
technique, such as 7 Cep and HD 41004, but none 
of them has a binary separation smaller than 20 
AU. The contrast between observations and theo¬ 
retical predictions implies that a close stellar com¬ 
panion not only decreases the planet formation ef¬ 
ficiency but also negatively influences planet evo¬ 
lution. iJang-Condelll (1201511 also discussed the ef¬ 
fect of mass ratio on planet formation. While we 
found the tentative correlation between hot Jupiter 
occurrence and equal-brightness stellar companion, 
it is difficult t o mak e a direct comparison because 
iJang-Condelll 1)201511 focused on stellar companions 
with separations smaller than 100 AU whereas all 
but one companions in our study have separations 


larger than 100 AU. 

8. SUMMARY 

We select a sample 84 KOls with 97 gas giant 
planets to study the influence of stellar companions 
on planet formation. We obtain AO images for 60 
KOls with two instruments: PHARO at Palomar 
200-inch telescope and N1RC2 at Keck 11 telescope. 
For the rest of the sample, we use AO images avail¬ 
able from the CFOP database. In total, we have 
detected 59 stellar companions around 40 KOls. 

Since some of the detected stellar companions are 
not physically associated with the KOls, we de¬ 
velop two methods of testing the physical associ¬ 
ation. The first method makes use of the color in¬ 
formation. For stellar companions with detections 
in multiple filters, we estimate their masses and ab¬ 
solute magnitudes based on their colors. Their dis¬ 
tances can be estimated from the absolute magni¬ 
tudes, apparent magnitudes, and extinctions. By 
comparing the distances of stellar companions and 
the KOls, we can test their physical association. 
The second method is based on the statistics of 
galactic stellar population. For stellar companions 
that are detected in only one filter, the color infor¬ 
mation is missing. Instead, we run galactic stellar 
population model to simulate the Kepler field. With 
the results, we estimate the relative probability of 
a gravitationally bound stellar companions and op¬ 
tical doubles/multiples. Then we can calculate the 
probability of physical association as a function of 
angular separation and differential magnitude of a 
stellar companion. With these two methods of test¬ 
ing physical association, we can effectively exclude 
the effect of foreground and background stars on 
stellar multiplicity rate. 

We find that the stellar multiplicity rate for planet 
host stars is 01 q% within 20 AU. In comparison, the 
stellar multiplicity rate is 18%±2% for the control 
sample, i.e., field stars in the solar neighborhood. 
The deficiency of stellar companions for planet host 
stars indicates that gas giant planet formation is 
suppressed by stellar companions within 20 AU. 
The stellar multiplicity rate for planet host star is 
34%±8% for separations between 20 and 200 AU, 
which is higher than the control sample at 12 %± 2 %. 
This suggests that stellar companions in this separa¬ 
tion range play a important role in gas giant planet 
migration. Beyond 200 AU, stellar multiplicity rates 
are comparable between planet host stars and the 
control sample. 

We explore whether stellar companions of differ¬ 
ent properties lead to different planet properties. 
We find evidence of three distinct populations of gas 
giant planets. They are separated by two charac¬ 
teristic orbital periods, 10 and 70 days. The stellar 
companions around stars with each planet popula¬ 
tion have different differential magnitude distribu¬ 
tions. For example, stars with HJs {P < 10 days) 
tend to have bright stellar companions with K band 
differential magnitudes smaller than 2 mag. Stars 
with gas giant planets with P > 70 days tend to 
have stellar companions whose differential magni¬ 
tudes peak at ^3.0 mag in K band, which cor- 
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respond to early type M dwarf companions. We 
emphasize that the three-population hypothesis is 
still tentative, and more follow-up observations are 
needed to either support or disprove it. 

If the three-population hypothesis survives future 
tests, it may have significant implications on gas gi¬ 
ant planet formation and evolution. First, the mi¬ 
gration of gas giant planets may be dependent upon 
the mass of a stellar companion. The majority of 
gas giant planets in our sample cannot form in situ 
because of a lack of building and accreting materi¬ 
als. They must migrate in via some mechanisms. 
The migration mechanisms could be the same for 
different planet populations, e.g., the Lidov-Kozai 
mechanism. Or the migration mechanisms could be 
different, e.g., the Lidov-Kozai mechanism for HJs 
and in-disk migration for longer-period planets. In 
either case, the migration mechanisms have to be 
consistent with the observed mass-dependency of 
stellar companions. 

Second, if HJs tend to be preferentially found 
in binary star systems with similar magnitudes, 
this population of HJs have been neglected be¬ 
cause of the concern of flux contamination. The 
HJs in WASP-94 A B (iNev eu-VanMalle et ^|2014f) 
and Kepler-13 AB ( Szabp et al.| 20111; Mazeh et ^ 


l2012t IShporer et al.l 20l4lJohnson et al.ll2014D svs- 

tems may be the tip of iceberg for this popula¬ 
tion. Surveys that do not bias against MSS can 
constrain the HJ occurrence rate in binary star sys¬ 
tems with similar magnitudes. Current and fu¬ 
ture space missions for transiti ng planets, e.g., the 
K2 r nission ( iHowell et al.ll201^ and the TESS mis¬ 
sion (|Ricker et al.ll20ll i7 may serve the purpose. In 
fact, the transiting signal of HJs in binary stars may 
already exist in ground-based transiting surveys. In 


addition, AO-fed spectrographs will play an impor¬ 
tant role in validating o r confirming HJs in binary 
stars (e.g., ICrennl[2nil . 
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Fig. 1. — The median 5 -(t contrast curve for AO images of 84 KOIs is shown in solid line. Dashed lines are 1-a deviation 
of the contrast curve. Detections within b" are shown as asterisks. When analyzing the detection completeness, each KOI is 
treated individually for the observation band in which the AO image was taken. A total of 59 visual companions around 40 
KOIs are detected (Table [IJ. 
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Fig. 2. — Typical completeness contours for 3 techniques used to detect and constrain stellar companions around planet host 
stars. The radial velocity (RV) technique is sensitive to stellar companions within ^^30 AU and with small or intermediate 
mutual inclinations to planet orbital planes (blue hatched region). Note that i is the angle between the stellar companion 
orbital plane and the sky plane, so z ~ 90° implies a small mutual inclination between the stellar companion orbital plane 
and a transiting planet orbital plane. Dynamical analysis (DA) is sensitive to stellar companions at larger mutual inclinations 
between the stellar companion orbital plane and a transiting planet orbital plane (green hatched region). The adaptive optics 
(AO) imaging technique is sensitive to stellar companions at wider orbits (red dotted region). The combination of these 3 
techniques contributes to a survey of stellar companions with high completeness. 




a [AU] 

Fig. 3.— Comparison of the cumulative stellar multiplicity rate between field stars in the solar neighborhood (blue) and 
gas giant planet host stars (red). Hatched regions represent 1-cr uncertainty regions. The stellar multiplicity rate for planet 
host stars is within 20 AU. In comparison, the stellar multiplicity rate is 18%di2% for the control sample. The stellar 

multiplicity rate for planet host star is 34%±8% for separations between 20 and 200 AU, which is higher than the control sample 
at 12%di2%. Beyond 200 AU, stellar multiplicity rates are comparable between planet host stars and the control sample 
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Fig. 4.— Results of K-S tests for the null hypothesis that stars with shorter-period and longer-period gas giant planets have 
stellar companions with similar distribution of differential magnitude. At each dividing period, we conduct a Monte Carlo 
simulation and show the median (solid), 68% (dashed) and 95% (dotted) confidence intervals. There may be three populations 
of gas giant planets whose formation and evolution are influenced differently by their stellar companions. The dividing periods 
for these populations are and ^^70 days. 



Fig. 5. — The distribution of differential magnitude for stellar companions. For stars with short-period gas giant planets 
{P < 10 days), stellar companions tend to be bright with small differential magnitudes. For stars with gas giant planets with 
P > 70 days, the differential magnitudes of their stellar companions tend to be fainter than the short-period counterparts. 
There is a peak at '^S.O mag, but it may be due to a lack of sensitivity for fainter stellar companions. 
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TABLE 1 

RV AND AO DATA FOR 97 KOIS 


KOI 

KIC 

a 

(deg) 

S 

(deg) 

KOI 

Kp Teff 

(mag) (K) 

logs 

(cgs) 

d 

(pc) 

#PL 

R.PL 

(R®) 

Period 

(day) 

^tart 

(MJD) 

RV 

T’end 

(MJD) 

KOOOOl.Ol 

11446443 

286.808472 

49.316399 

11.338 

5814.00 

4.380 

207.0 

1 

14.40±1.60 

2.470613 

54216.610352 

56208.22753 

K00003.01 

10748390 

297.709351 

48.080853 

9.174 

4766.00 

4.590 

38.8 

1 

4.68±0.18 

4.887800 

54336.352539 

56485.62695 

K00005.01 

8554498 

289.739716 

44.647419 

11.665 

5861.00 

4.190 

286.6 

2 

5.66±0.72 

4.780329 

54983.515625 

56486.43945 

KOOOIO.OI 

6922244 

281.288116 

42.451080 

13.563 

6025.00 

4.110 

919.7 

1 

15.90±2.10 

3.522499 

54983.540039 

55781.53418 

K00017.01 

10874614 

296.837250 

48.239944 

13.303 

5826.00 

4.420 

494.3 

1 

11.07±0.54 

3.234700 

54984.560547 

55043.51953 

K00020.01 

11804465 

286.243439 

50.040379 

13.438 

6011.00 

4.230 

608.4 

1 

17.60±2.50 

4.437963 

55014.412109 

55761.32519 

K00022.01 

9631995 

282.629669 

46.323360 

13.435 

5972.00 

4.410 

593.4 

1 

11.27±0.70 

7.891450 

55014.403320 

55792.43847 

K00046.01 

10905239 

283.255493 

48.355232 

13.770 

5764.00 

4.400 

604.1 

2 

4.33±0.37 

3.487688 



K00063.01 

11554435 

289.226166 

49.548199 

11.582 

5721.00 

4.470 

212.3 

1 

6.31±0.31 

9.434158 



K00094.02 

6462863 

297.333069 

41.891121 

12.205 

6184.00 

4.181 

463.6 

4 

4.22±0.42 

10.423707 



K00094.03 

6462863 

297.333069 

41.891121 

12.205 

6184.00 

4.181 

463.6 

4 

6.73±0.67 

54.319930 



K00094.01 

6462863 

297.333069 

41.891121 

12.205 

6184.00 

4.181 

463.6 

4 

11.40±1.10 

22.343001 



K00097.01 

5780885 

288.581512 

41.089809 

12.885 

5934.00 

4.040 

789.7 

1 

16.10±2.00 

4.885489 

55106.878906 

55115.87109 

K00098.01 

10264660 

287.708832 

47.333050 

12.128 

6395.00 

4.150 

446.2 

1 

10.00±1.40 

6.790123 

55440.500987 

56533.43359 

K00108.02 

4914423 

288.984558 

40.064529 

12.287 

5975.00 

4.330 

354.6 

2 

4.46±0.52 

179.601000 

55073.468750 

56145.49804 

K00119.01 

9471974 

294.559174 

46.062328 

12.654 

5632.00 

4.440 

313.0 

2 

3.90±2.60 

49.184310 



K00127.01 

8359498 

289.607971 

44.345421 

13.938 

5731.00 

4.450 

617.4 

1 

10.93±0.47 

3.578783 



K00128.01 

11359879 

296.200592 

49.140121 

13.758 

5786.00 

4.420 

579.8 

1 

11.97±0.85 

4.942783 

55284.481445 

55509.07519 

K00131.01 

7778437 

299.097534 

43.497589 

13.797 

6411.00 

4.400 

937.3 

1 

9.00±3.20 

5.014233 



K00135.01 

9818381 

285.240845 

46.668251 

13.958 

6082.00 

4.370 

810.1 

1 

10.56±0.46 

3.024095 

55752.008789 

55806.83398 

K00137.01 

8644288 

298.079437 

44.746319 

13.549 

5385.00 

4.430 

421.9 

3 

4.75±0.43 

7.641571 

55075.508789 

56146.48437 

K00137.02 

8644288 

298.079437 

44.746319 

13.549 

5385.00 

4.430 

421.9 

3 

6.01±0.54 

14.858940 

55075.508789 

56146.48437 

K00139.01 

8559644 

291.653168 

44.688271 

13.492 

5952.00 

4.380 

612.0 

2 

7.34±0.66 

224.797120 



K00141.01 

12105051 

288.038300 

50.651611 

13.687 

5425.00 

4.500 

463.1 

1 

5.43±0.29 

2.624234 



K00152.01 

8394721 

300.517120 

44.381580 

13.914 

6405.00 

4.420 

897.1 

4 

6.10±1.90 

52.091100 



K00157.03 

6541920 

297.115112 

41.909142 

13.709 

5685.00 

4.380 

514.8 

6 

4.18±0.76 

31.995467 

55440.500987 

56533.43359 

K00179.02 

9663113 

297.045410 

46.328701 

13.955 

6081.00 

4.420 

771.2 

2 

5.00±1.70 

572.397900 



K00244.01 

4349452 

286.638397 

39.487881 

10.734 

6103.00 

4.070 

313.5 

2 

6.51±0.89 

12.720365 

55366.602539 

56519.40820 

K00279.01 

12314973 

295.486511 

51.013500 

11.684 

6418.00 

4.280 

268.6 

3 

5.10±1.60 

28.454899 



K00289.02 

10386922 

282.945648 

47.574905 

12.747 

5812.00 

4.458 

348.4 

2 

5.04±0.00 

296.637100 

56449.400391 

56532.29101 

K00318.01 

8156120 

288.153992 

44.068821 

12.211 

6285.00 

4.290 

415.8 

1 

5.19±0.49 

38.583360 



K00319.01 

8684730 

290.117523 

44.872940 

12.711 

5952.00 

4.190 

421.3 

1 

7.50±1.10 

46.151590 



K00340.01 

10616571 

297.664673 

47.801392 

13.057 

5811.00 

4.400 

417.3 

1 

16.80±1.40 

23.673188 



K00344.01 

11015108 

283.340271 

48.549042 

13.400 

5984.00 

4.340 

597.5 

1 

3.80±1.60 

39.309240 



K00351.02 

11442793 

284.433502 

49.305161 

13.804 

6330.00 

4.430 

870.7 

6 

6.80±3.00 

210.596590 



K00351.01 

11442793 

284.433502 

49.305161 

13.804 

6330.00 

4.430 

870.7 

6 

9.80±4.20 

331.643000 



K00366.01 

3545478 

291.664185 

38.619255 

11.714 

6207.00 

4.180 

352.1 

1 

10.60±4.60 

75.112019 



K00367.01 

4815520 

284.472168 

39.911812 

11.105 

5569.00 

4.360 

153.6 

1 

4.98±0.71 

31.578680 



K00372.01 

6471021 

299.122437 

41.866760 

12.391 

5872.00 

4.487 

355.2 

1 

8.52±0.23 

125.630640 



K00375.01 

12356617 

291.201202 

51.144279 

13.293 

5757.00 

4.140 

778.2 

1 

10.40±1.40 

988.881118 



K00377.02 

3323887 

285.573975 

38.400902 

13.803 

5777.00 

4.450 

619.4 

3 

8.22±0.59 

38.907202 

55342.448242 

56506.36328 

K00377.01 

3323887 

285.573975 

38.400902 

13.803 

5777.00 

4.450 

619.4 

3 

8.28±0.59 

19.273938 

55342.448242 

56506.36328 

K00633.01 

4841374 

293.427032 

39.942429 

13.871 

6070.00 

4.030 

640.1 

1 

5.80±1.80 

161.479000 



K00638.01 

5113822 

295.559418 

40.236271 

13.595 

5980.00 

4.310 

575.9 

2 

3.80±1.20 

23.636883 



K00638.02 

5113822 

295.559418 

40.236271 

13.595 

5980.00 

4.310 

575.9 

2 

3.90±1.30 

67.093500 



K00672.02 

7115785 

291.169525 

42.640808 

13.998 

5524.00 

4.410 

524.6 

3 

4.03±0.41 

41.749990 



K00680.01 

7529266 

292.287323 

43.197281 

13.643 

6327.00 

4.350 

786.9 

1 

8.00±2.40 

8.600145 



K00682.01 

7619236 

295.197998 

43.269508 

13.916 

5592.00 

4.250 

632.7 

1 

10.80±1.40 

562.142400 



K00683.01 

7630229 

297.824310 

43.258381 

13.714 

5887.00 

4.390 

647.6 

1 

6.00±0.77 

278.122200 



K00686.01 

7906882 

296.840759 

43.647121 

13.579 

5559.00 

4.470 

494.5 

1 

11.40±4.40 

52.513565 



K00697.01 

8878187 

289.000885 

45.154270 

13.684 

5779.00 

4.050 

1315.7 

1 

4.24±0.75 

3.032154 



K00707.03 

9458613 

289.077545 

46.005219 

13.988 

5904.00 

4.030 

1236.7 

5 

4.10±0.36 

31.784660 



K00707.02 

9458613 

289.077545 

46.005219 

13.988 

5904.00 

4.030 

1236.7 

5 

4.69±0.41 

41.027950 



K00707.01 

9458613 

289.077545 

46.005219 

13.988 

5904.00 

4.030 

1236.7 

5 

5.50±0.48 

21.775754 



K00716.01 

9846348 

297.492157 

46.694519 

13.754 

6115.00 

4.490 

713.2 

1 

6.80±3.20 

26.893084 



K01162.01 

10528068 

288.868225 

47.759430 

12.783 

6138.00 

4.280 

490.7 

1 

4.40±2.30 

158.692400 



K01206.01 

3756801 

293.954590 

38.899971 

13.642 

5754.00 

4.110 

826.9 

1 

6.20±1.90 

422.917678 



K01274.01 

8800954 

283.256805 

45.087780 

13.354 

5310.00 

4.550 

356.9 

1 

4.73±0.28 

704.962626 



K01311.01 

10713616 

283.532959 

48.094261 

13.498 

6188.00 

4.200 

648.2 

1 

4.20±1.70 

83.577520 



K01335.01 

4155328 

291.020142 

39.220680 

13.968 

6222.00 

4.040 

1385.0 

1 

8.60±2.70 

127.832900 



K01353.02 

7303287 

297.465332 

42.882839 

13.956 

6260.00 

4.080 

1094.5 

2 

3.80±1.10 

34.543630 



K01353.01 

7303287 

297.465332 

42.882839 

13.956 

6260.00 

4.080 

1094.5 

2 

18.60±5.30 

125.865460 



K01375.01 

6766634 

288.320435 

42.261414 

13.709 

6169.00 

4.358 

594.0 

1 

6.78±0.00 

321.213900 



K01411.01 

9425139 

298.995087 

45.909512 

13.377 

5753.00 

4.410 

502.0 

1 

7.05±0.84 

305.056500 



K01431.01 

11075279 

287.022278 

48.681938 

13.460 

5649.00 

4.460 

486.8 

1 

8.45±0.48 

345.161300 

56472.390625 

56532.50195 

K01439.01 

11027624 

290.851776 

48.521339 

12.849 

5930.00 

4.090 

719.2 

1 

7.80±1.10 

394.610700 

55075.273438 

56531.31250 

K01474.01 

12365184 

295.417877 

51.184761 

13.005 

6293.00 

4.270 

552.0 

1 

9.30±1.20 

69.732970 



K01478.01 

12403119 

288.848633 

51.209049 

12.450 

5493.00 

4.417 

288.4 

1 

5.26±0.35 

76.133540 



K01645.01 

11045383 

298.210938 

48.559158 

13.418 

5197.00 

4.530 

379.9 

1 

10.31±3.10 

41.166759 



K01658.01 

4570949 

294.192108 

39.618999 

13.308 

6422.00 

4.320 

750.6 

1 

12.35±0.67 

1.544930 



K01684.01 

6048024 

293.534485 

41.329899 

12.849 

6387.00 

4.430 

616.1 

1 

7.30±2.60 

62.815570 



K01779.02 

9909735 

298.482819 

46.793621 

13.297 

5812.00 

4.140 

586.3 

2 

5.00±1.70 

11.815018 



K01779.01 

9909735 

298.482819 

46.793621 

13.297 

5812.00 

4.140 

586.3 

2 

5.80±1.90 

4.662723 



K01783.02 

10005758 

289.341431 

46.988239 

13.929 

6235.00 

4.460 

845.5 

2 

5.20±1.90 

284.042300 
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TABLE 2 

RV AND AO DATA FOR 97 KOIS 


KOI 

KIC 

a 

(deg) 

S 

(deg) 

Kp 

(mag) 

KOI 

T'eff 

(K) 

logs 

(cgs) 

d 

(po) 

#PL 

R-pl 

(R®) 

Period 

(day) 

^start 

(MJD) 

RV 

Tend #RV 

(MJD) 

K01783.01 

10005758 

289.341431 

46.988239 

13.929 

6235.00 

4.460 

845.5 

2 

8.50±3.10 

134.479720 



K01784.01 

10158418 

297.646057 

47.167488 

13.592 

5853.00 

4.540 

574.0 

1 

5.20±3.00 

5.007410 



K01792.01 

8552719 

288.971649 

44.624531 

12.160 

5689.00 

4.440 

300.1 

3 

4.61±0.27 

88.407030 



K01800.01 

11017901 

285.268585 

48.560009 

12.394 

5600.00 

4.430 

307.9 

1 

6.20±2.10 

7.794300 



K01805.02 

4644952 

288.811951 

39.770660 

13.828 

5708.00 

4.080 

789.5 

3 

4.30±1.20 

31.782260 



K01805.01 

4644952 

288.811951 

39.770660 

13.828 

5708.00 

4.080 

789.5 

3 

5.60±1.50 

6.941344 



K01808.01 

7761918 

294.743317 

43.461208 

12.487 

6278.00 

4.350 

424.4 

1 

4.00±1.60 

89.192840 



K01812.01 

6279974 

290.126526 

41.601082 

13.742 

6285.00 

4.420 

778.7 

1 

4.80±1.80 

0.805263 



K01825.01 

5375194 

295.300079 

40.556591 

13.895 

5545.00 

4.060 

904.3 

1 

4.20±1.10 

13.522604 



K02672.01 

11253827 

296.132812 

48.977402 

11.921 

5565.00 

4.330 

236.0 

2 

5.30±2.10 

88.516580 



K02674.01 

8022489 

289.651245 

43.824421 

13.349 

5973.00 

4.260 

501.1 

3 

7.30±2.60 

197.510340 



K02677.01 

9958387 

294.757690 

46.831120 

13.460 

6409.00 

4.340 

694.7 

1 

6.50±1.00 

237.788200 



K03444.02 

5384713 

297.429199 

40.561909 

13.693 

3842.00 

4.664 

122.4 

4 

5.74±3.20 

60.326632 



K03663.01 

12735740 

289.763611 

51.962601 

12.620 

6007.00 

4.340 

394.8 

1 

11.29±0.03 

282.525503 



K03678.01 

4150804 

289.791595 

39.285328 

12.888 

5650.00 

4.313 

413.8 

1 

9.12±0.04 

160.885644 



K03787.01 

7813039 

288.439117 

43.505249 

13.891 

5993.00 

4.355 

741.5 

1 

9.17±1.10 

141.733971 



K03791.02 

5437945 

288.474854 

40.651360 

13.771 

6340.00 

4.163 

932.2 

2 

5.65±0.21 

220.130023 



K03791.01 

5437945 

288.474854 

40.651360 

13.771 

6340.00 

4.163 

932.2 

2 

7.21±0.11 

440.785197 



K03811.01 

4638237 

286.153107 

39.714901 

13.906 

5551.00 

4.518 

505.5 

1 

6.98±0.99 

290.140253 



K03823.01 

4820550 

286.771454 

39.983822 

13.922 

5817.00 

4.544 

663.9 

1 

5.54±0.08 

202.117797 



K03875.01 

11911561 

290.149139 

50.239178 

13.579 

6022.00 

4.027 

660.4 

1 

4.73±0.70 

8.870420 



K03907.01 

7137213 

296.997803 

42.653320 

12.642 

6498.00 

4.081 

610.0 

1 

4.05±0.88 

28.643425 



K05515.01 

8429817 

291.452271 

44.431751 

13.952 

6247.00 

4.255 

806.9 

1 

8.88±1.50 

6.263490 
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TABLE 3 

Visual companion detections with AO data. 


KOI 

A Mag 

Separation 

Distance 

Detection 

PA 

Association 

ref. 




Primary 

Secondary 

Significance 


Probability 



(mag) 

(arcsec) 

(AU) 

(pc) 

(pc) 

('4') 

(deg) 




4.0 (i) 

1.13 



259.o14“1®2® 


135.0 


L14 


4.3 (r) 

1.11 





135.5 

0.98 

H14 

KOOOOl 

3.3 (z) 

1.11 


207.0122;4 



136.3 

0.99 

H14 

2.8 (J) 

1.12 

232.6 


233.5 

136.2 

1.00 

this work 


2.5 (H) 

1.11 

230.5 

207.0tlli 


383.6 

136.3 

1.00 

this work 


2.4 (K) 

1.12 

230.9 

207.0122-4 


68.5 

136.5 

1.00 

this work 

K00005 

2.3 (K) 

0.14 

40.3 

286.6111^ 


19.2 

308.9 

1.00 

CFOP 

KOOOlO* 

7.7 (J) 

3.04 

2795.9 

yxy. < i26.8 



94.3 

0.00 

A12 

KOOOlO 

6.2 (J) 

3.74 

3439.7 

yiy. /i26.8 


14.9 

89.3 

0.15 

A12 

K00017 

3.8 (J) 

4.01 

1982.1 

AQA 

4±y^.o g5 Q 


132.6 

39.9 

0.75 

A12 

K00020* 

7.9 (J) 

5.04 

3066.3 

608.41241“ 



139.6 

0.00 

A12 

K00097 

4.0 (J) 

1.90 

1500.4 

789.71??9“2 

3847.01^411.4 

49.8 

105.1 

0.90 

A12 


4.0 (K) 

1.91 

1508.3 

789.7lff9“2 


81.0 

105.1 

0.89 

A12 


0.1 (i) 

0.29 



707.41™ 


140.0 


L14 


0.5 (r) 

0.29 





140.0 

1.00 

H14 

K00098 

0.7 {z) 

0.29 





140.3 

1.00 

H14 


0.3 (J) 

0.27 

120.5 

446 

^^D.Z_2g Q 



143.7 

1.00 

A12 


0.4 (K) 

0.28 

124.9 

446.2ll°“„2 


15.7 

143.5 

1.00 

A12 

K00098 

6.2 (J) 

5.60 

2498.7 

446.2l4“12 

2604.514411 

20.8 

306.3 

0.13 

A12 


5.6 (K) 

5.59 

2494.3 

446.2ll°“o2 


12.0 

306.1 

0.16 

A12 

K00098 

7.2 (J) 

6.20 

2766.1 

446 

4^o.z_29 0 

422 ^ 

^zz.^_13g 4 

9.5 

237.9 

0.00 

this work 

6.3 (K) 

6.21 

2769.5 

446.2ll°“„2 


5.1 

237.9 

0.24 

this work 

K00108* 

7.2 (J) 

2.44 

865.2 

IsIelHy" 



74.9 

0.35 

A12 

K00108* 

7.2 (J) 

4.87 

1726.9 

354.61411 



112.4 

0.00 

A12 

K00119 

0.9 (i) 

1.05 



216.11^^14 


118.0 


L14 

0.2 (J) 

1.05 

327.5 

313.011^2® 


116.2 

118.4 

0.60 

this work 


0.2 (K) 

1.05 

327.5 

313.o11°®2« 


96.7 

118.4 

1.00 

this work 

K00137 

4.1 (J) 

5.59 

2359.2 

421 Q+4t).U 

4Zl.y 54 0 


116.7 

349.8 

0.57 

A12 

K00137* 

7.9 (J) 

4.80 

2025.1 

421 q+4t).u 

4Zi.y 54 Q 



340.5 

0.00 

A12 

K00137* 

7.5 (J) 

4.98 

2101.1 

421 



136.3 

0.00 

A12 


1.4 (i) 

1.10 



yo/.o_33g 4 


11.0 


L14 

K00141 

1.2 (J) 

1.06 

490.9 

463.11^11 


192.1 

13.9 

0.99 

A12 


1.4 (K) 

1.06 

490.9 

463.11^11 


242.4 

13.5 

0.99 

A12 

K00152 

5.7 (K) 

2.49 

2230.1 

897.li;“l“ 


6.0 

29.5 

0.23 

this work 

K00157 

4.4 (K) 

1.36 

698.7 

514.814414 


6.0 

179.9 

0.94 

this work 

K00157 

4.7 (K) 

4.09 

2104.3 

514.814414 


7.7 

238.0 

0.35 

this work 


3.5 (r) 

0.92 



457.01411" 


246.6 

0.99 

H14 

K00279 

3.1 iz) 

0.92 




247.2 

0.99 

H14 


2.3 (K) 

0.93 

250.4 

268.614®!® 


111.8 

246.9 

1.00 

CFOP 

K00340 

5.4 (K) 

5.34 

2227.7 

417 o+4y.» 


6.0 

57.6 

0.10 

this work 

K00344 

3.5 (K) 

4.13 

2470.2 

oy(.o i26.2 


45.9 

178.9 

0.68 

this work 

K00344 

5.2 (K) 

3.55 

2123.6 

cq7 c+^yu.u 

oy(.o 126,2 


9.7 

210.6 

0.34 

this work 

K00366 

6.5 (K) 

7.00 

2464.4 

352.i1®®7® 


5.6 

70.1 

0.02 

this work 

K00372* 

8.6 (K) 

2.49 

884.4 

355.21411 



157.8 

0.00 

A12 

K00372* 

8.0 (K) 

3.56 

1264.5 

ooo.z 54 Q 



56.9 

0.00 

A12 

K00372* 

8.2 (K) 

4.99 

1772.4 

ooo.z 54 Q 



170.7 

0.00 

A12 

K00372* 

4.0 (K) 

5.94 

2109.9 

355.21411 



32.7 

0.64 

A12 

K00375 

3.3 (K) 

5.47 

4254.0 

778.2l®|l5 


25.4 

157.0 

0.64 

this work 

K00375 

4.6 (K) 

3.19 

2486.2 

778.2l®|9l, 


5.9 

305.5 

0.56 

this work 

K00377 

4.5 (J) 

5.90 

3654.5 

(SI Q 4+^4-^ 

Diy.4_102,7 

2639.5l4'14 

37.5 

91.7 

0.31 

A12 


4.2 (K) 

5.89 

3648.3 

(SI Q 4+^'^-^ 
Diy.4 io2.7 


101.0 

91.7 

0.34 

A12 

K00377* 

6.8 (J) 

2.79 

1728.1 

(SI Q 4+^^-^ 
Diy.4_io2 7 

9717 

y(i/.o_j^745 g 


37.9 

0.10 

A12 

K00377 

6.6 (K) 

2.79 

1728.1 

619.4l®^2®7 


10.9 

37.8 

0.02 

A12 

K00633 

3.9 (K) 

0.67 

432.0 

640.ll®{l'‘ 


18.2 

18.4 

0.97 

this work 

K00683 

4.0 (K) 

3.42 

2214.4 

647.6l®14 


8.2 

268.9 

0.32 

this work 


Note. — *: stellar comoanions that are not detected bv our detection nioeline. References: A12 - lAdams et aT1 I2U1U): D14 -lUressine et al.l 120141') ; H14 
- IHorck et aTI ntTTll : L14 - ILaw et al.1 (Wm : LB14 - ILlllo-Box gTlH HUm) . 
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TABLE 4 

Visual companion detections with AO data. 


KOI 

A Mag 

Separation 

Distance 

Detection 

PA 

Association 

ref. 




Primary 

Secondary 

Significance 


Probability 



(mag) 

(arcsec) 

(AU) 

(pc) 

(pc) 

(c) 

(deg) 




0.0 (J) 

0.66 

868.4 

1 Qi c 7+209.8 

J-OJ-O. 1 707 Q 

f\7^ q + 1881.0 

0/i.O_424 2 

235.5 

55.3 

1.00 

this work 

K00697 

0.0 (H) 

0.66 

868.4 

-| 01 K 7 + 269.8 
1010 . < _727.9 


257.4 

55.2 

1.00 

this work 


0.0 (K) 

0.66 

868.4 

1 01 c 7 + 269.8 

lOiJ. t 727.9 


257.1 

55.7 

1.00 

this work 

K01274 

3.8 (i) 

1.10 


356.9titi 

412.21*1* 


241.0 


L14 


2.4 (K) 

1.09 

389.0 


9.0 

242.7 

0.99 

this work 

K01335 

4.6 (K) 

4.51 

6244.0 

1385.0+^;^«;^ 


17.3 

358.1 

0.26 

this work 

K01353 

5.4 (K) 

3.17 

3466.8 

1094.5+™^;“ 


7.0 

63.3 

0.19 

this work 

K01353 

5.9 (K) 

5.65 

6186.2 

1094.5+™^;“ 


5.0 

97.9 

0.00 

this work 


4.4 (i) 

0.77 


594.0^9^2® 

1907.7l“?l*-' 


269.0 


L14 

K01375 

3.8 (J) 

0.80 

473.2 


23.6 

270.0 

0.97 

this work 


3.6 (H) 

0.79 

467.5 



48.7 

269.5 

0.98 

this work 


3.6 (K) 

0.79 

470.3 

594.0tg4®2® 


26.4 

269.9 

0.98 

this work 

K01411 

5.3 (K) 

3.79 

1904.2 

502.0+^^;t 


7.1 

147.4 

0.14 

this work 

K01463 

6.4 (K) 

6.11 

2122.1 

347.5+73*4'’ 


5.4 

233.7 

0.03 

this work 

K01784 

0.9 (J) 

0.28 

160.7 

574.O+--7 

575.81“?!^ 

16.8 

288.4 

1.00 

this work 

0.8 (H) 

0.28 

160.7 

574.01?^+ 


10.6 

286.6 

1.00 

this work 


0.7 (K) 

0.28 

160.7 

574 


6.6 

291.0 

1.00 

this work 

K01808 

3.3 (K) 

4.68 

1986.0 



93.4 

163.1 

0.74 

this work 

K01812 

4.3 (i) 

2.37 


7170 7+168.5 
( ( 0 . (i 39.5 

2173.3l“ll 




LB14 


3.6 (K) 

2.37 

1844.2 


21.6 

297.9 

0.82 

this work 

K01825 

4.4 (K) 

5.56 

5025.4 

904.3+“*;* 


25.9 

295.3 

0.20 

this work 

K02672 

3.5 (K) 

0.69 

163.6 

236 . 0 +^^y 


44.2 

302.9 

0.99 

CFOP 

K02672 

5.9 (K) 

4.54 





308.0 


D14 

6.0 (K) 

4.61 

1088.0 

236 . 0 t^^® 5 '^ 


10.9 

310.4 

0.23 

this work 


2.8 (i) 

1.08 



179.4l?'ll 


9.6 


LB14 

K03444 

2.6 {z) 

1.08 


122 4+^^-® 
iZZ.q:_27 1 



9.6 

0.99 

LB 14 

2.2 (J) 

1.08 

132.2 


12.3 

9.5 

1.00 

this work 


2.4 (K) 

1.08 

132.2 

iZZ.4 27 1 


40.7 

9.5 

1.00 

this work 


4.5 (j) 

3.58 



1878.9l*??*e" 


264.4 


LB 14 

K03444 

4.7 {z) 

3.58 


122.4124-9 



264.4 

0.71 

LB 14 

5.0 (J) 

3.63 

443.8 


11.5 

264.8 

0.63 

this work 


5.3 (K) 

3.57 

436.7 

122.41219 


26.9 

264.8 

0.58 

this work 

K03678 

3.3 (K) 

2.61 

1081.7 

413 . 8 I 57 I 9 


33.9 

169.6 

0.85 

this work 

K03787 

5.1 (K) 

6.96 

5162.8 

741 


9.2 

254.9 

0.12 

this work 

K03823 

5.6 (K) 

2.33 

1546.0 

663.9lg**g* 


10.0 

58.0 

0.36 

this work 

K03823 

5.1 (K) 

5.06 

3357.5 

663.9l8**g* 


14.9 

239.4 

0.12 

this work 

K03907 

2.5 (J) 

2.77 

1689.5 

610.0ll'*g:2 

588.4l*“+2 

55.0 

74.2 

0.95 

this work 

2.2 (H) 

2.75 

1674.6 

610 

OiU.U_i35 7 


140.3 

74.3 

0.94 

this work 


2.1 (K) 

2.76 

1681.5 

diu.u^35 7 


97.2 

74.1 

0.96 

this work 


4.3 (J) 

1.59 

968.7 

610.0ll'fg;2 

QQ K + l^b/.d 

yy.o_29.o 

10.6 

162.3 

0.93 

this work 

K03907 

3.7 (H) 

1.57 

958.0 

eio.ol?!:? 

41.0 

163.4 

0.95 

this work 


3.5 (K) 

1.57 

955.2 

610.0l?l? 


27.9 

163.1 

0.96 

this work 

K05515 

4.1 (J) 

2.36 

1907.5 

806.91“*;? 

q+J-dODd.D 

zuoy.o_jg25 1 

24.3 

297.9 

0.79 

this work 

3.6 (H) 

2.36 

1907.4 

806.9lf3ll 


44.3 

297.7 

0.82 

this work 


3.7 (K) 

2.37 

1915.9 

806.9lf3l? 


23.2 

297.7 

0.80 

this work 

K05515 

5.4 (J) 

2.67 

2158.2 

806.9lf3|;? 

4443 

9.1 

114.2 

0.47 

this work 

5.1 (H) 

2.68 

2159.2 

806.9l?|l? 


12.9 

114.4 

0.45 

this work 


5.1 (K) 

2.67 

2155.6 

806.9l?«ll 


6.4 

114.5 

0.41 

this work 


Note. — *: stellar companions that are not detected by our detection pipeline. References: A12 - 1 Adams et al.| (2012|); D14 -|Dressins et al.| (2014|); H14 
-E orcl^^t_^^ {201^; L14 -[La^_et_a^ (201^ ; LB14 - [Lillo;Bo3^_et_a^ 







































